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We study the properties of classical and quantum compacton chains by means of extensive numer- 
ical simulations. Such chains are strongly nonlinear and their classical dynamics remains chaotic at 
arbitrarily low energies. We show that the collective excitations of classical chains are described by 
sound waves which decay rate scales algebraically with the wave number with a generic exponent 
value. The properties of the quantum chains are studied by the quantum Monte Carlo method and 
it is found that the low energy excitations are well described by effective phonon modes with the 
sound velocity dependent on an effective Planck constant. Our results show that at low energies the 
quantum effects lead to a suppression of chaos and drive the system to a quasi-integrable regime of 
effective phonon modes. 

PACS numbers: 05.45.-a,63.20.Ry,45.50.Jf 
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I. INTRODUCTION 

The investigation of nonlinear chains, started by Fermi, 
Pasta and Ulam in 1955 [l|, still remains an aetive and 
interesting area of research which attracts a significant 
interest of nonlinear community (see e.g. reviews |2t5| 
and Refs. therein). Usually, in such chains the nonlinear 
terms are relatively weak compared to the linear ones and 
strong nonlinear effects appear only at sufficiently high 
energy excitations. 

However, there are also other types of chains where the 
linear modes are absent and the dynamics is strongly non- 
linear at arbitrarily small energies. In such chains a time 
can be rescaled with energy and hence the system always 
remains in a strongly nonlinear regime. A prominent ex- 
ample is the Hertz lattice which describes elastically in- 
teracting hard balls where the elasticity parameter scales 
as a square-root of the displacement [6|-l9| corresponding 



to the nonlinearity index n — 5/2. Nesterenko |6l. llCHl2J 
described a compact traveling- wave solution in the Hertz 
lattice now known as compacton. A rigorous mathemati- 
cal description of compactons was given by Rosenau and 
Hyman [iSJ, |lj| for a class of nonlinear partial differential 
equations (PDEs) with nonlinear dispersion. A detailed 
analysis of compacton dynamics on lattices with various 
nonlinearity index has been performed recently in [15|. 
It was shown that on a finite lattice interactions between 
compactons lead to chaotic dynamics characterized by a 
spectrum of positive Lyapunov exponents. A well known 
example of such compacton chain is the toy "Newton's 
cradle" [H. 

In the compacton lattice the dynamics is chaotic at ar- 
bitrarily small energies when more than one compacton 
are exited. Therefore it is interesting to understand what 
happens in such quantum lattices and what are the prop- 
erties of quantum compacton vacuum. The studies of 



such quantum chains are reported in this work. We 
use the quantum Monte Carlo (QMC) method and the 
approach developed for the quantum Frenkel-Kontorova 
model as it is described in [l6[. The aim of this work 
is to understand the properties of quantum compacton 
vacuum and low energy excitations. It is interesting to 
note that the recent progress with cold atoms allowed 
to realize a quantum Newton's cradle \\j\ and to study 
energy redistribution between atoms. The experimental 
progress stimulated also theoretical studies of integrabil- 
ity and non-integrability in one-dimensional atomic lat- 
tices at high energy excitations [l8[- In contrast to that 
we study the properties of vacuum and low energy exci- 
tations in quantum compacton lattices when the linear 
terms are absent and nonlinearity is always strong in the 
classical case. Thus the classical dynamics of such com- 
pacton chain is always chaotic [15[ (except a special case 
of one compacton moving in a lattice). What are the 
properties of the quantum compacton chain is the sub- 
ject of studies of this work. 

The paper is organized as follows: the model descrip- 
tion is given in Section II, the properties of sound waves 
in the classical compacton chains are analyzed in Sec- 
tion HI, simple analytical estimates for the compacton 
chain are presented in Section IV, numerical results of 
the QMC are presented in Section V and the results are 
summarized in Section VI. 



II. MODEL DESCRIPTION 

The quantum compacton chain is described by the 
Hamiltonian 
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where index I marks M particles in the chain and n is 
the nonhnearity index. Here xi gives the particle co- 
ordinate counted from the equilibrium distance between 
particles which is taken to be a. For the Newton cradle 
a is given by the ball diameter. We use the dimension- 
less units in which the particle mass is equal to unity 
and the momentum pi gives the particle velocity. For 
the quantum problem the operators of momentum and 
coordinate have the usual commutator [pi, xi'] = —ihSi^i> 
with a dimensionless Planck constant h. We assume the 
periodic boundary conditions with xi = xi+m , Pi = Pi+M 
or fixed ends boundary conditions. The later is a specific 
case of even number of particles M = 2N. We note that 
the particles are distinguishable since they are located at 
well defined positions. 

Let us remind few known results for the harmonic chain 
at n = 2. By a canonical transformation to normal modes 

Hamiltonian ([Ij takes the form 
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At finite temperature T the relation ^ is modified for 
the usual expression for bosons 



h 

^k 



exp(fta;fe/T) — 1 



(5) 



With this form of ( S*^ ) one can obtain the expression for 

the correlator (xiXi+a) at finite temperature. 

The properties of the chain can be also characterized 
by the static form factor defined as 



F{q) 
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(6) 



where q can be viewed as a momentum transfer during a 
process of photon scattering, and a is the spacing of the 
chain lattice. 

Before to start the studies of the quantum compacton 
problem at n > 2 we consider in the next Section the 
propagation of sound waves in the classical chain at n ~ 
4. 



with the normal mode frequencies w^ — 2lo sa\{qk / 2) , 
u] ~ ^/a and wave numbers qt — irk/N, k — 1, . . . ,N~1. 
Here and below we use the fixed boundary conditions 
with xi=o = xi=N = and N particles. 

It is convenient to introduce sine modes via relations 
^k = Y :|- E; sin(gfcO*^ ^i = y "^ Efc sin(gfc/)S'fc. The 
vacuum state of the chain is a product of vacuum states of 
all modes. For any mode in the vacuum state one has an 



average of mode energy Uk being 
fvjjk/'^ and hence 



Uk 
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Thus different harmonics arc independent and as a 
result the squared deviation from equilibrium for 



a particle I is (xf) = jfJ2k=i sin^CftO 



-^ Efc=i sin^(gfc/)/a;fe. For the central particle / = N/2 
the displacement diverges logarithmically with the chain 
length (xf) = |rEEri'l/sin(qfc/2) « |ln(2/g^i„), 
where qmin = qk=i = i: /N . For h ~ 1 the displacement 
is of the order of unity for TV ^ 30. We will assume that 
such displacements are small compared to the distance a 
between particles. 

The spacial correlator of the chain can be also explic- 
itly calculated as (xiXi+a) = ^ J2k=i sin(gfeO sm{qkil + 
A)) (Slj where the brackets note the quantum average. 
Using ^ we obtain after summation over all I 



{xiXl+A)i 
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(4) 



III. SOUND WAVES IN CLASSICAL 
COMPACTON LATTICES 

To find the spectrum of sound waves and their de- 
cay rates we simulate numerically the dynamics of clas- 
sical system (jT]) at n = 4, a = 1 with the periodic 
boundary conditions. The simulations are done with a 
Runge-Kutta-Nystrom method with time step 0.02. Ini- 
tially a random distribution of momenta pi is seeded, so 
that the total energy per particle is 1. After a tran- 
sient of time Ttrans — 10 a perturbation of the type 
xi — >■ xi + s cos qkl is imposed with e — 0.1. The wave 
number qk is changed in the range < gfe < tt. The 
corresponding Fourier mode of xi changes in time ap- 



proximately as /„ ex e" 



"^1* cos ujqt. This time dependence 
is obtained from an ensemble of particles. We use up to 
10^ particles in an ensemble to obtain good averaging of 
statistical fluctuations. 

The numerical results are shown in Fig. [1] for the spec- 
trum of sound waves ui^q) and their decay rates 7(9). The 
frequency spectrum is close to the spectrum of sound in 
a harmonic lattice ui{q) = 2uismq/2 with ui w 1.4. At 
small q we have the spectrum of sound waves. The effec- 
tive velocity of sound for phonon like excitations is given 
by (Ii sa 1.4(a£';)^/'^. Here the dependence on a given av- 
erage energy per particle (Ei) appears since a typical fre- 
quency of particle oscillations is proportional to (Ei)^^^. 
The speed of sound in the physical lattice with distance 
a between particles is c = uJa. The decay rate of these 
waves drops algebraically with the decrease of the wave 
vector g as 7 ex g^. We obtain the value /3 « 1.69 that 
is close to the generic exponent /3 = 5/3 for the decay 
rate in nonlinear lattices (see e.g. [3[). Thus even if the 
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FIG. 1: (Color online) Bottom panels: spectrum of waves 
aj(g) in the classical compacton lattice ((T} at n = 4, a = 1 and 
energy per particle {Ei) equal to unity, at small wave numbers 
q the spectrum is close to a linear law of sound waves uj oc q. 
Top panels: decay rate of sound waves 7(5), at small q we 
have f (X q^ with /3 « 1.69. Left panels are in normal scale, 
right panels are in log-log scale, fits are shown by straight 
lines with indicated slopes. 



dynamics inside a compacton lattice is strongly chaotic 
(see [15[) the long wave oscillation properties of the whole 
lattice are well described by effective sound waves. In a 
certain sense the situation is similar to sound in a gas 
media: each particle moves chaotically but the collective 
long wave excitations are well described by sound waves. 
It is interesting to note that recently a localization of 
sound waves in a random three-dimensional elastic net- 
work of metallic balls has been observed experimentally 
in [l9| . Such a system can be viewed as a random three- 
dimensional Newton's cradle. However, our studies here 
are restricted to the one-dimensional case. 



IV. SIMPLE ESTIMATES FOR QUANTUM 
COMPACTON CHAINS 



The sound velocity c in a gas is given by the deriva- 
tive of pressure p over the gas mass density p at fixed 
entropy S (adiabatic process): c'^ = (g^) • Since the 

pressure is proportional to the force p oc dU/dxi, hence, 
c^ oc d^U/d'^xi. This leads to a simple estimate for 
the sound velocity based on the virial theorem accord- 
ing to which 2 (K) = n{U), where K and U are par- 
ticle kinetic and potential energies and brackets mark 
their average values. Since the temperature is propor- 
tional to the kinetic energy we have T ~ (K) = n{U) = 
a{{xi — Xj)") ^ a(Ax)" where Ax is an average dis- 
placement of a particle. This gives 



For the classical case this expression agrees with the 
above analytical results for w for n = 2, 4. 

For the quantum compacton chain we can use the 
Hcisenbcrg uncertainty relation p '^ H/Ax for the mini- 
mization of the ground state energy E = fi,^/2(Ax)^ -I- 
a(A.T)"/n that gives Ax - {h'^/a)'^^^''+'^\ and w^ - 
a(A.T)"-2 ^ Q,4/(«-2);j2(«-2)/(«+2)^ As a result the 

sound velocity of the quantum chain is 



C^ - a2^4/(Ti-2)^2(n-2)/(ri+2) 



(8) 



For n ~ 2 this velocity is independent of fi but for n > 2 
it decreases with H that corresponds to the decrease of the 
ground state energy of a nonlinear oscillator. The prop- 
erties of the ground state are analyzed in the numerical 
simulations presented in the next Section. 



NUMERICAL RESULTS OF QUANTUM 
MONTE CARLO 



A. Method description 

For our numerical simulations of quantum chain ([T]) we 
use the Metropolis algorithm (MA) [20[ in the Euclidean 
time T related to the system temperature T = fi/r. The 
simulations are done in the same way as for the studies of 
the quantum Frcnkel-Kontorova model described in [16| . 
A gen eral description of this QMC method can be found 
in^. 

The paths in the discrctizcd Euclidean time are gener- 
ated by the statistical sum 
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which links the problem to a statistical mechanics for 
configuration distribution of some lattice of TV x Nr size 
where N is number of particles and A^r = t j At is number 
of discrete steps of size At in the Euclidean time interval 
T. As usual the periodic boundary conditions are used 
in this time with xi(tj + t) = xi{tj) and tj = jAr. In 
our numerical studies we use up to Nr = 1000, r ^ 200 
and up to 2 • 10^ Metropolis updates. We fix a = 1 for 
numerical simulations. 

The numerical simulations assume a generation of con- 
figurations {xkj} with probability proportional to their 
weights in the statistical sum. The Metropolis method 
looks very efficient, providing an update step gives rather 
large modifications for a give site. However, correspond- 
ing modifications are local and are dominated by the 
nearest neighbor sites, that results in a significant slow- 
down for long wave configurations. Thus it is useful to 
combine the Metropolis method with the microcanonic 
dynamics (MCD) method. The MCD method is a noise- 
less algorithm and it works as follows: all variables xi,m 



arc considered as some coordinates, and the sum 
1 
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(10) 

as a potential energy of a certain system. Then, a set 
of auxihary momentum variables {Vi,j} is added and the 
equations of motion are solved numerically in an auxiliary 
update "time" variable u: 



&Pi,j/du 
dxij/du 



dU 

dxij 



Vi.-jUix) 



(11) 



With this dynamical description the system evolves over 
some iso-energy hypersurface in the phase space, and we 
get an ensemble of configurations. 
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2: (Color online) Typical data for the Fourier harmonics 
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FIG 

I ^9 , i^m I ^ E'lid its fit by the dependence (|14p ; black curve shows 
numerical data, dashed red/gray curve shows the fit (curves 
are overlapped). Data correspond to the linear chain at n = 2, 

ft = 1, Qj = Wtv/N, N = 64, Tmax = 100, Nr = 1000. 

However, an obvious disadvantage of such a method is 
that one needs to solve the differential equations numer- 
ically with a step which becomes smaller for decreasing 
At, due to the terms U{x) with Ar in the denomina- 
tor. But even worse, these terms act as some noise that 
reduces the relaxation rate along I (space dimension). 
Thus, in order to accelerate the relaxation processes we 
introduce into l/l(x) a parameter Ck 



U{x-Ck)^Y. 

1,3 



^^(x,,-x,,„i) +— (x, 



3 -Xl-X.]T 

(12) 
Then the update step is organized as follows: 

a) the fast mixing stage, with Ck = At, at which the 
smallness of denominator is canceled, here the MCD 
method is applied during the auxiliary "time" interval 
u - 100; 

b) return to Ck = 1 and application of the Metropolis 
updates at maximum 400 updates. 



Such a combined approach allows us to reduce signif- 
icantly the number of Metropolis updates and to have 
more rapid numerical simulations. We checked that both 
methods (only MA steps and MA steps combined with 
the MCD method) give the same results. 



B. Quantum excitations above the compacton 
vacuum 

The ensemble of quantum paths, obtained by the nu- 
merical methods described above, determines the proper- 
ties of the vacuum ground state of the system. Averages 
of various quantities over this ensemble give the corre- 
sponding average observables at this state. However, the 
study of fluctuations of quantum paths allows us to ex- 
tract a more interesting information about the spectrum 
of low energy elementary excitations. 

To extract this information we consider the Fourier 
harmonics of quantum paths 



Aj^'^™ = E "^^-'^ sin(gjZ) exp(ii^mT/c) 



(13) 



l.k 



where qj = ttj/TV, j = 1 
v,n = 27rm/Tmax, m = 0, . 



..,A^ - 1, Tfe = kAr, 
,Nt-/2. One expects a 
Lorentzian distribution in frequency to for exponential 
decay of quasiparticle excitations in the imaginary time: 
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Here 

2 



we use the renormalized frequency Dm = 
-^ sin(^'mAT/2) with the sine term appearing due to the 
discreetness of time steps. The fit of data for A, 



ij, 



al- 
lows to find the spectral dependence w = (jj{qj) and thus 
to determine the dispersion law of elementary quantum 
excitations. A typical example of such a fit is shown in 
Fig. [21 

The spectrum of low energy excitations extracted via 
such a procedure is shown in Fig. [2]for n = 2,4,8. For 
n = 2 the data reproduce the theoretical result for a 
harmonic chain with a) = uj(q)/(2sui{q/2)) — 1. For n = 
4, 8 this form of spectrum is preserved with a moderate 
renormalization of w » 1.2, 1.5 respectively. This shows 
that even if the classical compacton chain is fully chaotic 
the quantum compacton vacuum is rather regular and 
.is characterized by the phonon type excitations rather 
similar to the case of a harmonic chain. 

Our data show that w varies with h and n. This leads 
to the dependence of the sound velocity c = dio/dq on 
these two parameters. The numerically obtained depen- 
dence of c on ?i is shown in Fig. H) The fit by an alge- 
braic dependence c (x ff gives r] = 0.33; 0.56 for n = 4, 8 
respectively. These numerical values are close to the the- 
oretical power from Eq. ([5]) with 77 = (n — 2)/(n -|- 2) 
corresponding to 77 = 1/3, 3/5 for these n values. The 
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FIG. 4: The dependence of sound velocity c on h. Squares, 
triangles and circles correspond to n = 2, 4 and 8, re- 
spectively. The straight lines show the fit dependence (see 
text). Parameters of simulations are: A^ — 64, Tmax = 100, 
A^T- — 1000. Logarithms are natural. 
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FIG. 5: Spectra of quantum phonon modes for the quartic 
chain at n = 4. Data correspond to A'' = 64, r = 100, 
Nt = 1000, and to the interval h — 0.063 ~ 10, extrapo- 
lated smoothly from h = 0.063 to ft = 0. 



FIG. 3: (Color online) Spectra of quantum phonon modes 
shown by black bottom curve and obtained from QMC sim- 
ulations for n = 2 (top panel), n — i (middle panel) and 
n = 8 (bottom panel); top red/gray curve shows the ratio 
u}/ sm{q/2) for each panel. The data are obtained at ft = 1, 
N = 128, Nr = 1000, Tmax ~ 200, the simulations are done 
with MA using 2 ■ 10*^ updates. 



global variation of the whole spectrum Ld{q) with h is 
shown in Fig. [5] for n = 4. 

Trying to find deviations from a harmonic chain be- 
havior we compute numerically the amplitudes of phonon 
modes S^ (see Eq. ^ with q = k/{N — 1)) and the cor- 
relation function {xiXi+a)i (see Eq. (g])). The results are 
shown in Fig. [6l For n = 2 the numerical data are in 
a good agreement with the theory for a harmonic chain. 
For n = 4, 8 our numerical data show the dependencies 



rather similar to the case of a harmonic chain with slight 
vertical shift which can be attributed to the modified 
values of ui discussed above. We note that for the case 
of T > the theoretical formulas ([3|), ^ are computed 
with the expression ^ for bosons at finite temperature. 
Since the value of r in Fig. [S] is rather large there is no 
significant difference between the theoretical expressions 
for T = and T = H/t. 

An additional attempt to see deviations from a har- 
monic chain behavior is performed by computing the 
form factor F(q) of the chain given by Eq. ^. However, 
the results presented in the left panel of Fig. [7] show that 
all three chains with n = 2,4,8 give very similar curves 
for F{q) which practically overlap. 

The confirmation of similarity between these 
three types of chains in the ground state is 
given by the direct computation of the correlator 




FIG. 6: (Color online) Left column: Amplitudes of phonon 
modes, see Eq.©, for n = 2 (top), 4 (middle), 8 (bottom). 
Right column: Normalized correlation function, see Eq.Q, 
for I = 64 (central particle) for the same order of panels. 
Other parameters h = 1, N = 128, iVr = 1000, r = 200. 
Simulations include 2 • lO'^ updates. Red/gray and blue/black 
dotted curves give the theoretical expectations for harmonic 
chain with temperature T — h/r and T = 0, respectively 
(curves overlap). 



K{qi,q2) = (SgSg,) /{{SgSg) {Sg,Sg,)Y^^ which should 
be proportional to Sq^^q^ accordmg to ([3]). Indeed, the 
numerical data show a strong peak at qi ~ q2 with a 
residual noisy level of K at other qi ^ (72 without any 
structural dependence on qi,q2- This residual level can 
be characterized by the total weighted admixture of 
other modes to a given mode q via 



9' #9 



(15) 



This characteristic is shown in the right panel of Fig. [7] 
for n = 2,4,8. The admixture 'w{q) increases with n 
but still it remains rather small for strongly nonlinear 
lattices with n = 4, 8. This gives one more confirmation 
that the quantum compacton vacuum is rather close to 
a harmonic one. 

At finite temperatures the spectrum of excitations w(g) 
and the amplitudes of phonon modes S"^ are still well ap- 
proximated by the theoretical dependence for a harmonic 
chain as it is shown in Fig. [5] at T = 0.1 which is about 
two times larger than the excitation energy for a mode 
with minimal frequency huirnin = hir/N « 0.049. The 
situation is found to be qualitatively the same when the 
temperature is increased up to T = 1 keeping fixed other 
parameters of Fig. [5] even if at ji = 8 the splitting be- 
tween the theoretical curve and the numerical data be- 




FIG. 7: (Color online) Left panel: formfactor F{q), see Eq.®, 
for n = 2 (black curve), 4 (red/gray points), 8 (blue/black 
points) (data practically coincide). Other parameters are as 
in Fig. [6] including a = 1. Right panel: admixture to nor- 
mal modes from other harmonics. Red/gray and blue/black 
curves show data for quartic and octic chains, the dotted black 
curve corresponds to linear chain and gives an estimate for a 
noise level. Data correspond to h = 1, N = 64, r = 50, 
Nt — 500, number of independent quantum paths is 10^. 



comes more visible (due to a similarity of these data with 
Fig. [5] we do not show them here). 




FIG. 8: (Color online) Left column: spectra of quantum 
phonon modes at T = 0.1 for n = 2 (top) and 4 (bottom) 
with the curves as in Fig. [S] Right column: Amplitudes of 
phonon modes, see Eq.® in the same order; red/gray and 
blue/black dotted lines give theoretical expectations for har- 
monic chain with temperature T — h/r and T = 0, respec- 
tively. Other parameters are fi = 1, Nt = 200, Tmax = 10. 
Simulations are done with combined MA and MCD method 
with 10* Metropolis updates. 



VI. SUMMARY 

In this work we investigated the properties of collective 
modes in nonlinear compacton chains. For the classical 
compacton chains we find that the chain has sound modes 
which decay rather slowly due to nonlinear wave interac- 
tions with the rate 7 ex g*^ with (3 w 5/3. This decay rate 
is in agreement with the generic result of decay of sound 
waves in one-dimensional nonlinear chains |3|. 

On local scales the classical dynamics in such com- 



pacton chains is strongly chaotic. One could expect that 
this chaos may lead to nontrivial properties of quan- 
tized chains. However, our extensive numerical studies 
of quantum compacton vacuum show that it is character- 
ized by low energy phonon excitations which are rather 
similar to those of a harmonic chain. The main differ- 
ence is that the sound velocity of these phonon modes 
depends on an effective Planck constant as it is described 
byEq. dH). 

In a certain sense quantum effects suppress the signa- 
tures of classical chaos in the ground state. Such a phe- 
nomenon is known for quantum systems with a few de- 
grees of freedom. For example, a ground state of a Sinai 
billiard can be rather well approximated by a Hartree- 
Fock trial function with one maximum so that the signa- 
tures of quantum chaos appear only in the semiclassical 
regime for highly excited states |22| . This is more or less 
natural for systems with few degrees of freedom. Our 
case has infinite number of degrees of freedom but in spite 
of that the quantum compacton vacuum remains rather 
similar to a vacuum of a harmonic chain in which the 
oscillator frequency is dependent of an effective Planck 
constant. It is possible that certain signatures of such 
quasi-integrability at low energies find their manifesta- 
tions in a slow chaotization of excitations in a quantum 
Newton's cradle observed in experiments [l7| . However, 
we should note that such a statement can be considered 
only on a qualitative level of rigor since the model ([T|) at 
77. = 4 or 8 gives only an approximate description of ball 
interactions (see discussion at [15[)- 

However, let us note that our methods are well adapted 
for analysis of low energy oscillatory quantum waves. It 
is possible that other approaches should be used to de- 
tect quantum chock wave type excitations with large dis- 
placements between two parts of the chain (in principle 
the energy of such type compacton like excitations is not 
very high and is independent of the lattice size). Other 
methods should be used to analyze such type of excita- 
tions. 

It is possible that a degeneracy of chaos in a vicinity of 
quantum compacton vacuum is linked to the space homo- 
geneity of the model ((ij. Indeed, the ground state of the 
quantum Frenkel-Kontorova model has much more rich 
properties [16| appearing due to a presence of periodic 
potential. 
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